Wave communication across regular lattices 
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We propose a novel way to communicate signals in the form of waves across ad- dimensional 
lattice. The mechanism is based on quantum search algorithms and makes it possible to both search 
for marked positions in a regular grid and to communicate between two (or more) points on the 
lattice. Remarkably, neither the sender nor the receiver needs to know the position of each other 
despite the fact that the signal is only exchanged between the contributing parties. This is an 
example of using wave interference as a resource by controlling localisation phenomena effectively. 
Possible experimental realisations will be discussed. 

PACS numbers: 03.67.Hk, 03.65Sq, 03.67.Ac, 42.50.Ex 



Localization phenomena in linear wave systems are 
closely linked to wave interference effects. Anderson lo- 
calisation in disordered media is a prime example thereof 
still posing challenges to both theory and experiment 
[l[ 60 years after its discovery @. Recently, a new re- 
search area has emerged focusing on interference as a 
resource and making use of localisation phenomena in a 
controlled way. Prominent examples are among others 
time reversal imaging [H and reconstructing the Green 
function in terms of correlation functions 0], see also @ 
. Here, information about the wave system is obtained 
by manipulating a seemingly 'noisy' signal using phase 
coherence. We will focus here on another class of wave 
localisation phenomena with counterintuitive properties, 
namely (quantum) search algorithms and (quantum) ran- 
dom walks. Wave search algorithms gained prominence 
with Grover's work Q demonstrating a \J N speed-up 
compared to a classical search within an unsorted data 
base of N items. Even though search algorithms became 
an inherent part of quantum information theory Q, the 
speed-up is in effect caused by wave interference as has 
already been pointed out by Grover @] and has been im- 
plemented for a classical wave system in Q. Based on 
ideas from quantum random walks [IoI - H^ ]. Grover's al- 
gorithm has been generalized to spatial search algorithms 
on networks such as on a hypercube [13 . [H[ and on reg- 
ular lattices [l5| . Experimental realisations of quantum 
random walks have been achieved again both using clas- 
sical waves (optics) [l6[ and quantum devices [l7| . 

Starting from wave search algorithms, we will 
demonstrate that localisation can be used to establish 
communication channels across a regular lattice with 
surprising properties: (i) signals can be exchanged 
exclusively between a source and a receiver point, where 
neither the sender nor the receiver know the position of 
each other; (ii) the signal can track a moving receiver 
in the network; (iii) the algorithm can be used as a 
searching device without the necessity to know the time 
of measurement, (a typical requirement for Grover's 
search algorithm); (iv) the protocol can act as a sensitive 
switching device for wave transport through a lattice; 
(v) the algorithm can be effectively implemented both 
on a quantum computer and using classical waves 
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FIG. 1: Regular grid with d — 2, n = 8; local scattering within 
a unit cell at vertex x is described by the matrix a. 



only. We will first describe the set-up of the search 
algorithm. We then introduce a simplified model for 
the search and explain the wave communication protocol. 

We consider wave propagation across d-dimensional 
periodic lattices of identical scatterers or periodic po- 
tentials with fixed lattice parameter a, see Fig. [TJ It is 
important that the lattice has a finite number of sites n 
along each axis with a total number n d of lattice sites. 
To simplify the calculations we will restrict ourselves to 
models with nearest-neighbour interaction only and con- 
sider periodic boundary conditions. The wave dynamics 
within each unit cell is given by a local scattering matrix 
a mapping incoming channels onto outgoing channels, 
sec Fig. [T] (a is also denoted a coin matrix in the con- 
text of quantum walks). The overall wave dynamics is 
then given in terms of an operator Uq mapping incom- 
ing onto outgoing wave coefficients between unit cells. 



We have dim U 



where m denotes the number of 



open scattering channels within a unit cell. Furthermore, 
Uq is unitary when disregarding dissipation. Stationary 
solutions are obtained by the condition 



det (1 



_ifca 



U o )=0 



(1) 



where k is the wave length and exp(ifca) is a phase 
shift between incoming and outgoing waves. We ne- 
glect any (in general weak) k dependence of a and 
thus Uq. Note that the spectrum obtained from ([1]) is 
now periodic in k with period 27r/a; the eigenvalues are 
kjj = (2wl — 0j)/a, I S N , where 6j are the eigenphases 
of Uq with j = 1, . . . , N. To start with, we will consider a 
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model consisting of a single open channel between nearby 
lattice sites (m = 2d) and we assume Kirchhoff bound- 
ary conditions at each vertex. The Hilbert space is then 
effectively N = 2dn d dimensional. This physical model 
captures the essence behind the effect described below. 

We label incoming wave components from each unit 
cell as 1^) <£> \x) = 1^,5;), where x specifies the vertex 
in position space and i , for i = 1, . . . ,d gives the ± 
direction in dimension i. Incoming waves at vertex x 
are mapped onto outgoing waves by a scattering matrix 
a. For Kirchhoff boundary conditions, one obtains 
a = 2\s) (s| — lad an d I s ) is the uniform distribution 
I s ) = ^Ei=i + I*")) 0- Outgoing waves in 
direction are now identified with incoming waves 
at an adjacent vertex x ± ae.i where efj is the unit 
vector in direction i. The local scattering processes 
is described in terms of a (global) scattering (or coin) 
matrix C = a <8> l n <». The full wave propagator (or 
quantum walk) Uq is obtained from C after identifying 
incoming and outgoing waves of adjacent unit cells 
accordingly. The spectrum of the unperturbed walk Uq 
exhibits a band structure as shown in Fig. [5] a) , here for 
d = 2. (For a finite lattice, the quasi-momenta k x , k v are 
discretised according to K x ,y = 27rj/ (na); j = . . . n — 1.) 

Following Ambainis, Kempe and Rivosh (AKR), the 
quantum walk Uq acts as a search algorithm after mark- 
ing a target vertex \v) by a modified scattering matrix 
a' [l5| . that is, one considers C = C — (a — a') ® \v) (v\. 
The AKR search uses a' = —lad- Since \s) is an eigen- 
vector of a, we may write U 1 = Uq (1 — 2 \sv) (sv\), where 
\sv) = \s) <£>\v). The search algorithm is initialised in 
the uniform state |$o) = V \sx), and the walk 
(U') T |$ ) localizes at |u) after T cx vOV steps. 

The AKR search can be analysed by defining a one 
parameter family of unitary operators 



Ux = U + (e I7rA - 1) Uq \sv) (sv\ 



(2) 



one obtains Uq for A = or 2 and the AKR search for 
A = f . The part of the eigcnfrcqucncy spectrum of U\ in- 
teracting with the perturbation is shown in Fig.[5]b). The 
spectrum is periodic in k with period 27r/a independent 
of A. When varying A, a "perturber state" \v\) emerges 
which crosses the k mod 27r/a = axis at A = 1. The re- 
sulting avoided crossing between the initial state |$ ) an d 
\v\) is shown in Fig.[5]c). Note that |$o) is the fully sym- 
metric eigenstate of Uq corresponding to a d-dimensional 
Bloch-vector k — of the unperturbed spectrum with 
eigenvalue k mod 2ir/a = 0. Like in Grover's algorithm, 
the quantum search U' = U\=i rotates the initial state 
|$o) into a localised state \v\) which has here a strong 
overlap with the target state \sv). The search time T is 
inversely proportional to the gap at the avoided crossing 
A, that is, T w n/A. 

In order to obtain an estimate for the search time To 
as well as the efficiency of the search, that is, the matrix 
element (sv\v\), it is essential to find the approximately 




FIG. 2: a) The band structure at A = for d = 2 and in 
the limit n — ¥ oo with k, the Bloch wave numbers; b) the 
eigenphases of U\ for n = ll,d=2;c) the avoided crossing 
with spectral gap A at A = 1 and k mod 2n/a = together 
with the approximate eigenstates \i>\) and |$o) (dashed lines). 



invariant two-level subspace near the crossing s pan ned 
by |3>o) and |i/a=i). The technique developed in [14|] for 
the hypercube has been adapted to regular grids. We will 
only give the result here, further details will be presented 
elsewhere uM. One finds for the normalised vector |^a=i) 



K=i) = -<H«a=i>\/^X>-^x (3) 

K#0 



where v is the position of the target vertex and 
and are the eigenvectors and eigenphases of the un- 
perturbed walk Uq [Ta ] . The d-dimensional label n with 
Ki = 0, . . . , n — 1 is equivalent to the (discretised) Bloch 
wave number, see Fig. [5] a). The eigenphases are ex- 
plicitly given as cos#g = 3X4=1 cos ^if 1 - The overlap- 
matrix element (sv\i/\ = i) can be estimated 181 1: 



0{l/yf\ogN) for d = 2, 
0(1) for d > 2. 



Detailed expressions for the leading order coefficients for 
d = 2 and 3 are given in [18|. We find that | i^a= 1 } is ex- 
ponentially localised on the marked vertex v; note, that 
the overlap of \sv) with a typical eigenstate of the unper- 
turbed spectrum is of the order 0(N~2) <C 0(1). 

Near the avoided crossing, the level dynamics can be 
described in terms of the two-level sub-space spanned 
by the orthogonal vectors \v\) and |$o)- Writing the 
unitary operator U\ as U\ = e~ lHx , one obtains at A = 1 
in the {|$o), |^x=i)} basis an effective two-dimensional 
Hamiltonian H\ = i of the form 
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T = 



T = 15 




T = 30 



T = 45 




FIG. 3: Probability distribution of the quantum walk on a 
31 x 31-grid up to T = 45 time steps. 



with ki = 2i:l/a,l £ N being the wave numbers corre- 
sponding to k = states of the unperturbed lattice and 
e is a real and positive coupling parameter, that is, 



-rMw«^+or) ( 5) 



with A, the gap at the avoided crossing. 

The start vector (1,0) = |$q) is rotated into the lo- 
calised state (0, 1) = |f.x=i) m T) = 7r/(2e) steps leading 
to the y/N speed-up [ll|. The whole process is 2 To- 
periodic , that is, one needs - like for Grover's algorithm 
[f| - to know the period To to perform the search. For a 
simulation of the search on a 31 x 31 grid, see Fig. [3J 

Interesting applications emerge when considering sev- 
eral target vertices, ,i = 1, . . . m with m <C N. We 
now define a set of parameters A = (Ai, . . . , A m ) and a 
search algorithm of the form 



U: 



-it/ 



(6) 



At A = (1, 1, . . . , 1) and ka mod 2ir = 0, one finds that 
there are m — 1 degenerate eigenvalues and two further 
eigenvalues forming an avoided crossing with the degener- 
ate subset. The corresponding set of m+1 eigenstates co- 
incides in good approximation with the subspace spanned 
by the uniform distribution |$q) an d now m localised 



states i/jA ,i = 1, 



Each of the |^) is well de- 



scribed by the approximation (|3| and {v\\v 3 \) ~ The 
localised states interact at the crossing predominantly via 
|$o) which takes on the role of a carrier state. In analogy 
to ((3]), we write a model Hamiltonian at the crossing in 



the basis {|$o) 



'A=i/ 



A: 



=i» as 



jj(m+l) X (m+1) 
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Like for the full propagator U\ = i, the spectrum of 
2j(m+i)x(m+i) cons i s ^ s f ( m _ eigenvalues equal to 

kia and two eigenvalues kia ± \pme with eigenvectors 
|cj ± ) = i/V2m (=Fi\/m, 1, . . . , 1) . The gap between the 
non-degenerate levels is now A = 2y / me with e given in 
([S]). Starting the search in the totally symmetric state 
|$o) at A = (1,...,1), one finds that £/J|<I>o) localises 
on all m marked vertices after To ~ n/2\/N/m (or, for 
d = 2, T ~ tt /2y / N/mlogN) steps simultaneously. 

More interestingly, the quantum walk can also be used 
to transmit signals across the network. Such a sender- 
receiver configuration has to the best of our knowledge 
not been described before and may have interesting appli- 
cations both in a quantum setting, but also for classical 
waves (such as microwaves, in optics or acoustics). In- 
stead of starting the walk in the uniform state \$>o), we 
propose to begin the walk at one of the localised states 
|^™), say. At the avoided crossing, this state can approx- 
imately be described as 



K) 



+ V2(m-l)|wd)) (8) 



with | wo) = (m(m — l)) _1 / 2 (0, 1 . . . , 1 — m)* in the basis 
spanning i{( m + 1 ) x ( m + 1 ) ■ | Wo ) i s a vector in the degener- 
ate eigenspace with eigenvalue kia. Applying the walk 
for T s = n/(c\/m) = 2T steps leads to 
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This implies that a signal of intensity 4/m 2 is transmit- 
ted from the sender (located at the m-th marked ver- 
tex) to each of the m — 1 other marked vertices. The 
intensity at the sender at time T s is then of the order 
(1 — 2/m) 2 . These findings have been verified numerically 
in our model. Interestingly, in the case m = 2 with a sin- 
gle receiver, the signal is transmitted in full . This opens 
up the possibility of transferring signals directly between 
two points on a network where neither the sender nor 
the receiver know each other's position. In addition, the 
sender has information about the number of receivers by 
recording the signal at time T s . 

The effect persists also for a continuous source at the 
sender position. In Fig. 01 we recorded the signal both at 
the sender and at the receiver. (To keep the signal finite, 



4 




time 



FIG. 4: Sending continuously: signal at sender (full) and 
receivers (dashed/dotted). The response of the system to 
switching the position of the receiver at T = 200 is shown. 
(We introduced a damping of 1% at all vertices.) 

we added a small amount of absorption across the net- 
work). Again, without prior knowledge of the receiver's 
position, the network localises at the two marked vertices, 
thus making it possible to actually exchange information 
continuously between these two points. Changing the po- 
sition of the receiver leads to a sudden drop of the signal 
at the old receiver position and a build-up at the new po- 
sition, see Fig. 01 The system is thus capable of tracking 
a moving receiver position! The signal speed is limited 
by the transfer time T s ~ -JN (T s ~ y/N In N for d=2) 
and thus is the speed at which the receiver can move. 

The continuous sender/receiver protocol can also be 
used to search a marked item without knowing the search 
time To; this is a serious complication of Grover-type 
search algorithms when the precise number of marked 
items is not know (as Tq depends on to). Here, we find 



the marked items as long as we wait for times T > T s . 
Furthermore, the system can act as a switching device. 
Wave transport between two points on the grid can 
only be achieved, if the system is tuned to the avoided 
crossing. Slight detuning by for example changing the 
parameter A in ([2]) will quickly cut-off the signal. 

The described effects open up completely new ways 
of transmitting signals across regular networks. While 
a construction of the map U is certainly feasible in a 
quantum setting and can be implemented efficiently on a 
quantum computer 0] , an implementation using classical 
waves may be even more promising. For dispersionless 
wave dynamics and on regular lattices, the unitary 
matrix U\ is equivalent to a discretised version of the 
time dependent Green function (for times t — a/c with 
c being the wave velocity; corresponding time scales in 
a quantum setting would be given by the group velocity 
[201]). Indeed, localised states due to a local perturbation 
in a regular lattice are well known; (for example in 
optical crystals, see [IH). We predict that the effect 
will occur if "defect states" created by local phase 
perturbations (equivalent to the perturbed coin C) are 
pushed into the (discretised continuum) of the band 
close to the fully periodic state - the k = state, see Fig. 
[5] a. Furthermore, the interaction between the defect 
states and the lattice states must be small enough to 
lead to avoided crossings between these two states only. 
We expect that a signal (such as a laser or a microwave 
transmitter coupled into a periodic structure at a defect 
position) can be transmitted and focused onto another 
defect in the same lattice using the described effect. The 
totally symmetric state |$o) acts then as a carrier state 
guiding the signal between the perturbations. 
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